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ABSTRACT 

By means of a multimass isotropic and spherical model including self-consistently 
a central intermediate-mass black hole (IMBH), the influence of this object on 
the morphological and physical properties of globular clusters is investigated 
in this paper. Confirming recent numerical studies, it is found that a cluster 
(with mass M) hosting an IMBH (with mass M,) shows, outside the black hole 
gravitational influence region, a core-like profile resembling a King profile with 
concentration c < 2, though with a slightly steeper behaviour in the core region. 
In particular, the core logarithmic slope is s < 0.25 for reasonably low IMBH 
masses (M, < 10 _2 M), while c decreases monotonically with M m . Completely 
power-law density profiles (similar to, e.g., that of collapsed clusters) are admit- 
ted only in the presence of a black hole with an unrealistic M. > M. The mass 
range estimate 12s — 4.8 < log (M./M) < —1.1c — 0.69, depending on the mor- 
phological parameters, is deduced considering a wide grid of models. Applying 
this estimate to a set of 39 globular clusters, it is found that NGC 2808, NGC 
6388, M80, M13, M62, M54 and Gl (in M31) probably host an IMBH. For them, 
the scaling laws M. ~ 0.09(M/ M Q ) a7 M and M. ~ 50(cr obs /km s -1 ) 1 - 2 M , 
are identified from weighted least-squares fit. An important result of this 'collec- 
tive' study is that a strong correlation exists between the presence of an extreme 
blue horizontal branch (HB) and the presence of an IMBH. In particular, the 
presence of a central IMBH in M13 and NGC 6388 could explain why these clus- 
ters possess extreme HB stars, contrarily to the 'second parameter counterparts' 
M3 and 47 Tuc. 

Key words: black hole physics - stellar dynamics - methods: analytical - 
methods: numerical - globular clusters: general - stars: horizontal branch 



1 INTRODUCTION 

The presence and origin of intermediate-mass black 
holes (IMBHs) at the centre of globular clusters, are 
still being debated inside the astr ophysical community 
(see the recent g eneral reviews in Ivan der Marell |2004| ; 
iRasio et al.l l2006h . There are essentially two main for- 
mation theories: they could be Population III stel- 
lar remnants (see, e.g., iMadau fc Reed bOQll ) or could 
form in a runaway mergin g of young stars in sum 



ciently dense clusters (e.g.. iPortegies Zwart et al.l 12004 



Giirkan. Freitag. fc Rasidl2004l ; lFreitag. Giirkan fc Rasiol 
20063). 



The IMBHs existence is suggested primarily by an 
argument of plausibility: they would fill the 'embarrass- 
ing' gap between stellar black holes (BHs; with mass 
M m < 10 Mq) and super-massive BHs (having M, ~ 
1O 5 -1O 9 M ) that are believed to reside in the nucleus 
of most galaxies (e.g.. lKormendvil2004l : lRichstonell2004l : 



F 



err arose 



fc Fordl 
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for recent reviews). By exploit- 
ing the accurate kinematical measurements of gas and 
stars in nearby galaxies, important correlations between 
the mass of super-massive BHs and various global prop- 
erties of their hosting galaxy have been deduced (see 
iFerrarese fc Fordl [20051 . and references therein). For in- 
stance, extrapolating the lMagorrian et alj (|l99fet) relation 
- M, ~ 10 -3 times the mass of the host system - to glob- 
ular clusters, these would contain BHs with an intermedi- 
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ate mass M, ~ 10 2 -10 4 . Indirect observational evidences 
come from the detection of ultra-luminous X-ray sources 
radiating at super- Eddington luminosity (> 10 39 erg s _1 ) 
and thought to originate from matter infall on BHs con- 
sider ably more massive than stellar BHs (see iFabbianol 
120061 . for a review on this subject). 

In globular clusters, however, direct kinematic obser- 
vations are seriously obstacled by crowding and by the 
relatively small number of stars inside the BH gravita- 
tional influence region (hereafter BHIR), which cannot 
be well resolved in many cases. Thus, so far, only very 
few clusters have been the object of accurate kinematic 
studies so as to infer about the possible presence of mas- 
sive objects in the central regions. 

We just remind t he cases o f the core- 
collapsed clusters M15 llGerssen et all 12001 120031 ; 
iMcNamara. Harrison fc Anderson! 120031 ) and 47 Tuc, 
in which such a presenc e is not confirmed yet by the 
latest observational data llvan den Bosch et al ] |200rj and 
iMcLaughlin et all 120061 . respectively), and the highly 
concentrated and massive cluster Gl (in M31), where a 
central ~ 2 x 10 4 Mp IMBH has b een claimed to reside 
l|Gebhardt. Rich, fc Hdliool . l2005h . while a recent deep 
photometric analysis of the core region of u) Cen seems 
to be consistent with an IMBH with mass ~ 10 4 Mq 
l)Novola. Gebhardt fc Bergmannl 120061 ) . Note, however, 
that for both M15 and Gl, direct and accurate iV-body 
simulations reproduce the observed mass-to-light ra- 
tio wi thout the need of an IMBH ijBaumgardt et aD 
l2003al lbl), although the peculiar scenario of a two 
clusters merging event is required in the case of Gl. 

Other, less direct, insights can be gained from the 
time behaviour of the period of millisecond pulsars sited 
in clusters central region, which permits to deduce the 
local acceleration. This occurs in the intriguing case of 
NGC 6752 where there are indications of the existence 
of an unde rluminous and compact component w ith mass 
~ 10 3 M l|Ferraro et al.ll2003l : IColpi et al.ll2005l . and ref- 
erences therein). 

It is interesting to note that the above-mentioned 
individual studies concerned globular clusters that are 
supposed to host IMBH because of their (present) high 
central density and (presently) high rate of stellar colli- 
sions in their compact cores. Nevertheless, present condi- 
tions may be drastically different from those at the early 
ep ochs of clusters life. More imp o rtantl y, as firstly argued 
by iBaumgardt. Makino. fc HutJ (I2005T ). their dynamical 
status could be even inconsistent with the long-term ef- 
fects that a IMBH actually produces on the cluster in- 
ternal evolution. Using collisional iV-body simulations, 
IBaumgardt et all found that the high stellar density in 
the vicinity of the central IMBH enhances the rate of close 
encounters that, in turn, induces a rapid expansion of 
the central region, giving rise to a medium-concentration, 
King-like profile, whose features are almost independent 
of the BH mass. They argued that an observable finger- 
print of the IMBH presence is just a slight slope of the 
density in the core region, thus suggesting as probable 
candidate clusters hosting IMBH, those having a surface 
brightness logarithmic slope ~ —0.2 in the core (while the 
typical projected profile of core-collapsed systems goes 
liker" 1 ). 



Similar conclu sions have been recently reached by 
iTrenti et al.l l|2007h who conducted direct TV-body exper- 
iments of clusters with IMBHs and realistic fractions of 
primordial binaries: the expansion of the core is further 
enhanced and the formation of the density cusp is con- 
fined in the immediate vicinity of the BH. This leads to 
density profiles with a core to half-mass radius ra tio sig- 
nifica ntly higher than when the IMBH is absent l|Trentil 
2006). All these studies are telling us that a IMBH repre- 
sents a 'heat' source in the central region, acting in a way 
similar to harden ing binaries in contra sting gravothermal 
collapse (see also lHeggie fc Hutll2003l ). 

However, besides these global morphological hints, 
other fingerprints of the IMBH presence could be 
revealed. As known, one of the most important effect 
of the gravitational influence of the IMBH on the 
surrounding stars, is the tidal erosion and disrup- 



Frank & Reesl 19761; Amaro-Scoane, Freitag & Spurzem 


20041; IBaumgardt. Makino. & Ebisuzakil 


2004a; 


Freitag, Amaro-Scoane & Kalogera 


2006a; 



Baumgardt et al.l 120061 ). Besides complete disaggre- 
gation, this may lead to a loss of envelope mass of 
passing-by giants. The mass loss from stellar outer 
layers is also one of the possible explanation for the not 
yet ascertained origin of 'blue subdwa rfs' - also called 
extreme horizontal branch stars, see iRich et all Il997l 
and, for a review on this subject in t he context of dense 
galactic nuclei, also lAlexa nder 2005], section 3.4. These 
are stars heavier and bluer than turn-off stars, and they 
are located at the lower-left end of the horizontal branch 
(HB), in the region of the color- magnitude diagram 
referred to as extreme HB (hereaf ter EHB; see the 
discussion in iMevlan fc Heggidll997l . sect. 9.7, and the 
recent review in lCa tclan 2003). 

Given the cuspy behaviour of the stellar density in 
the BHIR and the relatively high ratio between M, and 
the single stellar mass, it is reasonable to expect that the 
above-mentioned mechanism of blue subdwarfs formation 
is significantly enhanced by the presence of the IMBH. 
Thus, the possibility to find some connection between 
this presence and that of an abundant population of EHB 
stars naturally emerges and, to our knowledge, has never 
been explored so far. This is one of the aim of the present 
paper. 

While the direct investigation of the dynamical evo- 
lution of a cluster harbouring a massive compact object 
is still at the beginning - because of the recent devel- 
opment of sufficiently advanced computational tools - 
the study of the stellar phase-space distribution function 
(DF) around a massive BH in the quasi-steady regime is 
a classical subject in Stella r Dynamics tackled by several 
authors since the early 70s dPeebleJ 19721; iBahcall fc Wold 



Sh 

Icoi 



1976 ; IShapiro fc Lightman 19761; iLightman fc Shapiro 



1977i ; lCohn fc Kulsrudlll978T ). If a time much longer than 
the cluster relaxation time is passed since the BH forma- 
tion and accretion, dynamical and 'thermal' equilibrium 
can take place and a Maxwellian DF ~ exp(— E/a 2 ), 
with E the total energy and a a velocity parame- 
ter, should be established. Nevertheless, close enough 
to the compact object, stars are either disrupted, be- 
cause of the strong tidal interaction, or swallowed by 
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the BH, thus generating an outward flux of positive en- 
ergy that, in stationary conditions, must be constant 
and uniform through the central region. This precludes 
the Maxwellian distribution b eing valid in the BH vicin- 
ity. Simple scaling arguments ijBinnev fc Tremainelll987l . 
section 8. 4.7) or treatments b ased on the 'heat-transfer' 
equation (|Freitag et a l. 2006b, section 2.2) show, indeed, 
that a realistic space density behaviour is ~ r~ 7 ^. 

Such a power-law density (or 'cusp') corresponds 
to a DF ~ (— -E) 1//4 , as was first found, through 
ap proximate soluti o ns of the Fokke r-Planck equation, 
by iBahcall fc Woll 1 19761) (hereafter IBW76I ) and more 
recently confirmed with the help of more accurate 
methods incl uding relaxation effects, such as: M onte- 
Carlo codes (|Shapirol Il985l; iFreitag fc Bend [2002]) . gas- 
dynamical methods ( Amaro-Seoane et al.l [2004) and 
direct iV-body simulations (|Baumgardt et al ] l2004al : 
IPreto. Merritt fc SpurzemlbOoi T 

However, equilibrium self-consistent models of stel- 
lar systems including massi ve BHs are no t so nu mer- 
ous. We just mention the iGerssen et al.l (|2002l ) pa- 
per on Ml 5 employin g spherical non-parametric mod- 
els and the pape r s by Ivan den Bosch et all l|2006h and 
Ivan de Ven et ah! (|2006l ). in which axisymmetric models 
of, respectively, M15 and uj Cen are generated to infer the 
radial behaviou r of the M/L ratio, b y means of a gener- 
alization of the ISchwarzschildl (1 19791 ) orbit superposition 
method. A further generalizatio n of this technique has 
been also recently presented in ICapuzzo Dolcetta et al.l 
(2007) to model triaxial galaxies with dark matter and 
a central density c usp. It is worth noting, finally, that 
iMcLauehlin et all (|2006l ) . in order to derive hints about 
the presence of an IMBH in 47 Tuc, applied a single- 
mass and isotropic King model, even if in a not fully 
self-consistent way. 

In the present study, a self-consistent parametric 
model is constructed and used to estimate the mass of 
possible IMBHs in clusters and their influence on the 
cluster density profile, as well as on the presence of EHB 
stars. It basically consists of a multimas s King model 
appropriately extended to incorporate the IBW76I DF in- 
side the BHIR. The model is spherical, isotropic (in the 
present version) and relatively easy to construct and use. 
Its multimass nature makes straightforward the inclusion 
of central energy equipartition, which is suggested to play 
a non-negligible role in the presence of a massive BH, too 
(IBahcall fc Wolfj Il977l; iMurphv. Cohn, fc Durisenl Il991 



Baumgardt. Makino. fc Ebisuzakil " 2004bl ; Freitag et all 
2006bl; lHopman fc Alexanderll2006h . 



The model will be described in section [2] while in 
section [3] the observable features of the generated sur- 
face brightness and velocity dispersion profiles will be 
illustrated, also with respect to the debated structural 
properties of clusters hosting IMBH. Then, in section [4] 
the model will be employed to deduce the mass of possi- 
ble IMBH s from the central de n sity p rofiles prov ided by 
the recent iNovola fc GebhardtJ (|2006l ) (hereafter iNGOfj ) 
analysis of WFPC2 Hubble Space Telescope (HST) obser- 
vations of 38 galactic globulars (Gl in M31 is considered 
as well). In section [5j finally, the connection between the 
existence of IMBHs and the presence of EHBs will be 
investigated. Conclusions will be drawn in section [6] 



2 DESCRIPTION OF THE PARAMETRIC 
MODEL 

In the following, given the length-, velocity- and density- 
scale parameters, r , a and p, respectively, it is convenient 
to deal with dimensionless quantities, namely: 

radius x = r/r; 
velocity w = v/o~; 
mean total potential W (x) 
energy (per unit mass) E \ 
density v(x) = p(x)/p; 

cluster mass in stars inside the radius x, M(x) 
M(x)/pr 3 . 

BH mass p = M,/pf s . 



= -<f(x)/a 2 ; 
w 2 /2-W; 



A suitable DF can be given by generalizing the King 

(|l966l ) lowered Maxwellian, so as to include the IBW76I 
DF below a proper 'transition' energy, i.e. 

f{E) = { (2^)- 3 /V 
0, 



1). 



if E < -Wbh, 

if -W BH ^ E < 0, (1) 

if E > 0, 



where the coefficient c = (2-k) ' (e 



3/2/„W BH 



1)W, 



-1/4 



sures the / continuity and W is defined so to have 
W(x t ) = where x t is the model 'limiting radius'. The 
'transition' potential Wbh = W(xbh) is considered to be 
the potential on the surface of the sphere, with radius 
ibh, coinciding with the BHIR. 

The BHIR can be defined as the largest sphere, cen- 
tered at the BH position, in which the motion of the stars 
is dominated by the field generated by the compact ob- 
ject. This definition is coherent with the assumption of 
Maxwellian behaviour just outside this region (i.e., for 
E ^ —Wbh). Then, a reasonable condition for xbh is 
that the enclosed mass in stars is a small fraction of the 
BH mass, i.e. 



-M(:ebh) = O.l^i, 



(2) 



even if this forces to solve (as we will see later) an 
implicit equation that involves a mass profile that is 
only a-posterior i known . Thi s condition is similar to 
that adopted in iMerrittl (|2004l ): nevertheless we chose a 
smaller fraction (0.1 instead of 2) to ens ure the validity 
of the power- law behaviour for f{E) fsee IBW76I . section 
IIIc). 

Th e usually adopted formu la ibh ~ GM./(v 2 }f 
(|BW76l : iBaumgardt et ail l2004al lbh is not a-priori com- 
patible with this BHIR definition, because it gives a ra- 
dius that yields typically M(xbh) < W~ 4 p. Indeed, even 
if the BHIR has to be sufficiently small for the BH grav- 
itational field to dominate the dynamics inside, on the 
other hand it must be large enough to make plausible 
the hypothesis of isothermal distribution outside it. This 
latter condition would not be realistic in an environment 
in which t he BH still exerts a strong influence on the 
dynamics (|BW76l ). 

The DF of equation ^ guarantees the 'dynamical 
equilibrium' of the system, in the sense that, according 
to the Jeans theorem, it is a solution of the collisionless 
Boltzmann equation. The presence of a massive BH, as- 
sumed to be located at rest in the cluster centre, makes 
that the potential time-independence is preserved and 
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thus the DF of equation |T]) is still a valid equilibrium so- 
lution, because E remains an integral of motion. There- 
fore, one can proceed in the same manner as in obtain- 
ing usual King models, that is - as required by the self- 
consistency -just by solving, for W(x), the dimensionless 
Poisson's equation in spherical symmetry 



d 2 W 2 dW 

dx 2 x dx 



-(4TrGpr 2 a~ 2 )u(W), for x > 0, (3) 



where the origin is excluded to avoid the BH singularity 
and the stellar density, expressed as a function of W , is 



v{W) 

where 
ui(W) ~ 
and 

u 2 {W) = 

Vl(W B 



= 4ty 



f(E)w 2 dw = 



za, if W < Wbh, 
u 2 , if W > Wbh, 



(4) 



e^erf (</W) 



4tyc 



W 



" 2 + l W 3/2 



(5) 



W 



w 
~2 



1/4 



w dw + 



V 2W BH 

2 7/2 TTcg(W)W 7/4 + za(Wbh) 



(6) 



Choosing xo = ibh, these conditions define a two- 
parameters (Wbh, Wbh) set of modelifl with all the oth- 
ers being just scaling parameters. 

By the definition of BHIR it can be assumed that 



. , GM. 

dr (xBnr) 2 



x 



— P r > 

BH 



that, using equation ([9]), leads to 



4nxl H ' 



(12) 



(13) 



Thus, the whole set of models can be conveniently de- 
scribed by the pair (p, "Wbh), while the solution inside 
the BHIR can be found by integrating equation <[3j back- 
wards from ibh to a given minimum radius, which we 
generally set equal to Q.Ixbh- 

Finally, ibh must be a solution of equation ([2} with 
M (x) given self-consistently by the model that uses ibh 
itself as a BHIR radius. This is found iteratively, evaluat- 
ing the 'new' radius Xb™ through equation ([2| with M (x) 
generated employing the previous value of ibh- The pro- 
cess starts with an initial guessed ibh = GM,/3a 2 r 
and is halted when the difference between two subse- 
quent values stays within a given tolerance (we chose 



Nbh -xbh\/x bh < 0.1). 



The first term of the r.h.s. of equation ((6} is the density 
corresp onding to a polytropic st ellar model with index 
= 7/4 |Binnev fc Tremainelll987h . with 



g(w) 



I ■ 3/2 

Ism ' y 



7/2 , 



') dy, 



(7) 



where 6 = cos 

The integral in equation (0 yields 



g(w) 



4 
21 



F 



F 



2.1 Inclusion of a mass spectrum 

The continuum stellar mass spectrum of the real cluster 
can be represented by a set of n stellar components, each 
one having a star mass mk, a total mass and a velocity 
parameter at- The DF of the more realistic multimass 
case, /mm, can be given as a linear combination of the 
f of equation (fp (se e, e.g., |Pa Costa fc Freeman! 1 19761 ; 



iGunn fc Grimrj|l97' 



f mm (E) = J2 <*kfk{E) = J2 VkfWkE), 



(14) 



— (5 cos 6 + 3 cos 30) J sin 9 = 
42 



4 
21 



'S-§)-'(5)+^-^-«) V4 ^ 



,(8) 



where ui = W BH /W and F((f>) = f*(l - 2sm 2 y)-^ 2 dy 
is an elliptic integral of the first kind that can be accu- 
rately evaluate d by standard nu merical iterative proce- 
dures (see, e.g.. |Press et al.|[l988l . sect. 6.11). 

As in standard King models, it is convenient to im- 
pose the relation 



_ 2 = 9a 2 
r 4nGp 

With this choice equation (0 becomes 

d 2 W 2dW n .„.. . 

—r^- H 7— = -9v(W), for x > 0, 

dx 2 x dx 



(9) 



(10) 



that has a family of solutions depending only on the 
boundary conditions (with xo > 0), 



W{x ) = W . 



(11) 



where otk serve to reproduce the given mass function and 
Pk=a 2 /a 2 . 

A more rigorous treatment would require to 
explicitly depend on irik through the exponent of E 
in the DF inside the BHIR, as a consequence of 
mass segregatio n . Thi s dependence was quantified by 
iBahcall fc WolJ (|l977f ) and then confir med (at least 
parti ally) by m eans of Fokker-P l anck llMurphv et all 
Il99ll ), iV-bodv (iBaumgardt et all 1200401 ) and Monte- 
Carlo (|Freitag et alj 1200681 ) simulations. Nevertheless, 
the small gain in 'theoretical coherence' might be vanified 
by the use of a discrete set of stellar masses instead of the 
real continuous mass spectrum. Furthermore, the model 
would be complicated by the need of separate numerical 
quadratures in equation (f3]) for each mass class. For these 
reasons, we simplified the model construction by consid- 
ering the same 1/4 energy exponent in f(E < —Wbh), re- 
gardless of the stellar mass. However, the energy equipar- 
tition, responsible of mass segregation, can be still repro- 
duced at the border of the BHIR by a suitable choice of 
Pk (see below). 



In standard King models Wq = at xq = 
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Thus, the multicomponent Poisson equation has the 
same form of equation (|10p apart from v that is replaced 
by the total density 



I/tot(W) 



a 3/2 
fc=1 Pk 



(15) 



where v is given here by equation (j4|. 

Once a component is arbitrarily chosen as the 'refer- 
ence' one, say the first one, without loss of generality one 
can assume a = ai, i.e. f3i = 1. Then, since the set of ak 
are constrained by the given set of Mk, the solutions of 
the Poisson equation are determined just by the bound- 
ary conditions at a;o = ibh- Indeed, the set of /3k>i 
can be fixed by the requirement of energy equipartition in 
the cent r al reg ion. To do that, the procedure described in 
iMiocchil (|200rJ ) is applied at the radiufl xbh- In practice, 



once given a and VKbh, we have 
'mi 



Pk = W^k 



BH, 



(16) 



with k the function defined in equation (8) of iMiocchl 
j20rah . and o 2 k = o 2 /p k . 

Thus, Wbh and fi determine a 2-parameters family 
of models [through equation (|13[l ] also in the multimass 
case. Note, furthermore, that does now depend on m k , 
through /3k- 

As regards ibh, the same iterative procedure de- 
scribed in the single-component case is employed to find 
the value that solves equation ([2]). Finally, as a consis- 
tency test, the velocity dispersion profile generated by 
the model is compared with t hat deduced by the spher i- 
cal isotropic Jeans' equation l|Binnev fc Tremaind 09871 ). 



(■w')j(x) 



27 



4w t ot (a;) 



iyt t(y){M(y)+fi]y- 2 dy, (17) 



and the model output [M{x),Vtot(x)\ is accepted only if 
the difference between the two profiles (evaluated as a \ 2 
summation over the spatial grid points) is small enough 
(e.g., to within 90 percent of confidence level). 



3 SURFACE BRIGHTNESS PROFILES 
WITH IMBH 

To understand how the presence of an IMBH influences 
the shape of a typical projected profile, it is instructive to 
examine the results of a single- mass model first. In Fig.sQ] 
and [2] various projected surface density profiles (E) are 
plotted, fixing Wbh = 9.5, for some of the models written 
on Table [T] As a comparison, a standard King model is 
also plotted with Wo = 9.5. 

As expected, we can see that the BH gives rise to a 
cuspy behaviour extended from the centre up to a radius 
that can be called the 'cusp' radius, x cu , above which the 
much flatter core region begins. An 'operative definition' 
of this radius can be given as the radius where the surface 
brightness profile shows a concavity change (while x cu = 
if the inflection point is not evident, i.e. if the cusp is 
not observed). However, this radius was found to be well 



Because only for x > £bh all f k are lowered Maxwellians. 



Table 1. Relevant parameters for single-mass models gener- 
ated with Wbh = 9-5. Some of them are plotted in Fig.s [T] 
and[2] fi is the dimensionless BH mass, M./M is the BH mass 
in units of total cluster mass, c is the concentration defined 
according to Equation il8l , £'bh i s the BHIR radius, x cu is 
the cusp radius, x c is the core radius and s is the logarithmic 
slope of the projected surface density around x cu - In the last 
line Xcu > xt (i.e. c < 0) indicates the presence of a full cuspy 
profile (in this case x c is a meaningless quantity). 
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0.78 


0.14 


14 


66 


7.7 x 10" 


-2 


1.8 


1.8 x 10" 


2 


0.77 


0.14 


84 


330 


0.12 




2 


9.2 




< 


0.14 


> x t 




> 1 





approximated - at least in the range of BH masses studied 
here (M./M < 2 per cent) - by the radius containing half 
the BH mass in stars [M(x C u) = At/2]. 

Moreover, it is convenient to define the 'core' radius 
x c as equal to the radius at which the surface brightness 
drops to 1/2 the value at x C u- Indeed, it is important to 
adopt a concentration parameter that is unambiguously 
defined when a density cusp is present. In particular, the 
parameter 



c = lo; 



(£) 



(18) 



is preferable to the standard King concentration param- 
eter ck = log(a;t), because this latter is based on the 
length-scale r that has no immediate relation with any 
relevant morphological feature of our IMBH modelqj. 
However, it has been checked that an IMBH model with 
a given c yields a profile that is well matched, in the 
core and in the tidal region, by that of a standard King 
model with ck = c. Note, also, that in the profiles de- 
picted in Fig.s [T] and [21 x c corresponds to the location of 
the 'knee' of the curves, in analogy to a standard King 
profile. Thus, in clusters whose central luminosity cusp 
is not resolved, one can infer about the presence of an 
IMBH directly using the standard ck parameter coming 
from observations. 

From these profiles (plotted as a function of the ra- 
dius measured in unit of x c ) three relevant features can 
be immediately noticed: 

(i) the projected density in the core region is not as 
flat as in standard King model, (see inset in Fig.s [TJd and 
®>) 

(ii) the presence of a central IMBH reduces signifi- 
cantly the concentration of the profile in comparison with 
a BH-free model with Wo = Wbh, and 



3 On the contrary, in King models x c ~ 1, i.e. f ~ the core 
radius. 
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x/x c 

Figure 1. (a) Behaviour of the potential and (b) projected 
surface density profile (scaled so as to give £ = 1 at the cusp 
radius) for single-mass models with Wbh = 9-5 and with var- 
ious IMBH mass as indicated (/i = means a standard King 
model, with Wo = 9.5). The dots mark the BHIR radii, while 
the arrows indicate the cusps radii (x cu ). The panel (b) inset 
shows the profiles in the core region around x cu . The concen- 
tration parameters are also reported in panel (a). 

(iii) the concentration decreases for increasing fi. 

These features are confirmed (as we will see later) in the 
case of multimass models, too. 

The last two points can be easily understood from 
the behaviour of the potential (Fig.s and that, 
at ibh, has a different slope when compared with the 
(practically flat) King model. The potential with the BH 
decreases below Wbh, then becomes flatter in the core 
(before dropping in the tidal region) leading the system 
to a lower concentration profile in comparison with the 
BH-free case. This suggests that the presence of a IMBH 
induces profiles with a large core. 

Note, from Table [1] the rapid growth of the cusp 
radius for /i > 1; for /i = 2, it becomes even larger than 
the limiting radius making, in fact, £(ic) to "jump" to a 
completely cuspy behaviour (not shown in the Figures), 
without any recognizable core. Since this leads to M, > 
M, such kind of configurations are discarded throughout 
our analysis, as unrealistic for globular clusters. Therefore 

- once the validity of the DF of equation JTJ is accepted 

- reasonable values of the BH mass at the centre of GCs 
are compatible only with a core-like profile. The density 
cusp is confined in the very inner region of the core. 

All this is in remarkab le ag reement with the re sults 
of lBaumgardt et all (|2005h and iTrenti et al l l|2007l ). ob- 




0.01 0.1 1 10 

x/x c 

Figure 2. The same as in Fig. [T] for models with higher BH 
mass. Since in this case the cusp radii are nearly equal, only 
one arrow is plotted. 



tained by means of collisional iV-body simulations in- 
cluding IMBH. These simulations have shown that the 
presence of an IMBH initially enhances the rate of ex- 
change of energy in close encounters among stars moving 
in the cusp region; this gives rise to a strong cluster ex- 
pansion which, in turn, yields final profiles well fitted by 
King models, with a rela tively low concentration (see also 
iBaumgardt et alj|2004al lbh. 

Finally, the apparent contradiction that can be seen 
in Fig.s HH2] of a decreasing BHIR radius for increasing 
BH masses, is actually due to the scaling of the radius 
with x c and to the fact that the core radius grows quite 
rapidly while ibh is nearly constant (see Table [lj . The 
decrease of the ratio ibh/ic is a consequence of the lower 
and lower concentration that these models show for grow- 
ing BH masses. Another consequence is the unexpected 
decrease of the logarithmic slope (s, with E ~ of the 
surface density around a; cu , that can be seen for fj, > 0.5 
(see Table [T]and Fig. [2b inset). 



3.1 The multimass case 

At this point, one could study in more detail the depen- 
dence on the BH mass of both c and the slope of £(ie) 
in the core region. Nevertheless, it is more appropriate 
to make this analysis directly in the multimass context, 
which is (for globular clusters) a much more realistic sce- 
nario than that of the single- mass case, because the mass- 
segregation could heavily change the effect of the pres- 
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Figure 3. (a) Concentration parameter obtained for a grid of 
values of M./M, and VKbh; no realistic models can be gener- 
ated above the dashed line, (b) Behaviour of c as a function 
of BH mass for some values of WbH; as indicated. 



ence of the BH in the brightness profile ([Baumgar dt et all 
120051 1. 

Thus, a stellar mass spectrum is included assum- 
ing a Salpeter mass function (diV oc m _1 ' 35 d log m) and 
considering 7 mass bins in the range 0.25 - 0.83 Mq, 
along with l ight and heavy re mnants, following the pre- 
scriptions of lCote et alj (|l995t ). Nevertheless, their light- 
est mass class was not included because it is below the 
lo wer mass limit predicted, under energy equipartition, 
in iMiocchil (|2006h . See Table [2] for the list of parameters 
for the adopted stellar components. 

The qualitative behaviours of the apparent concen- 
tration and of the £ slope as a function of the BH mass 
discussed in the single-component case, are basically con- 
firmed by the curves plotted in Fig.s[3p-[4p. They refer to 
a grid of models computed in the range 15 < Wbh ^ 9.5 
and 10 -2 ^ (i ^ 0.15. In this case the concentration has 
been evaluated on the total projected surface brightness 
profile, I(x), and the logarithmic slope s is again such 
that / ~ x~ s . It is worth noting that: 

(i) c decreases monotonically with the BH mass 

(ii) according to our model, IMBHs with M. > 10~ 3 M 
can exist only for non-collapsed distributions (c < 2) ex- 
hibiting a logarithmic slope in the core which is < 0.25 
for reasonably low BH masses (M. < 10 _2 M) 

(iii) s increases for growing BH mass up to a maxi- 
mum value, then the lower concentration achieved makes 
the core region to be much wider than the cusp region, 
causing the decrease of s. 

These features fully confirm the iBaumgardt et ail (|2005f ) 
findings. In particular, the logarithmic slope is compati- 
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Figure 4. As in Fig. [3]for the logarithmic slope of the surface 
brightness at the cusp radius. 



ble with the range of values found for the clusters studied 
by these authors' iV-body simulations. 

In Fig. |3p it is apparent a certain degree of self- 
similarity of the various curves; of course, lower concen- 
trations can be in principle obtained starting from lower 
Wbh - though this imposes an even higher mass cut-off 
of the stellar mass function at the low-mass end in or- 
der to ensure energy equipartition. Thus, no lower limit 
on c can be fixed. On the contrary, an upper limit to 
the concentration exists (roughly indicated by the dashed 
line in Fig. The same argument holds for s (Fig. 2}, 
for which lower values could be obtained for lower Wbh- 
Therefore, combining the upper limits of the two quan- 
tities, a range of admitted values is found for the BH 
mass: 

11.6s - 4.85 < log (^) < -1.14c - 0.694. (19) 



3.1.1 The effect of mass segregation 

The presence of mass segregation (as due to the imposed 
energy equipartition at sbh) is evident in Fig. [5] Three 
components are represented: the brightest component, 
the heavy remnants one and the lightest MS stars com- 
ponent (k = 1,2 and 7, respectively, in Table [2}. The 
first component gives almost all the cluster luminosity. 
The model has been generated with M^bh = 9.5 and 
\i = 0.155, corresponding to M./M = 6 x 10~ 3 . 

In the core region the lightest (heaviest) stars 
have the highest (lowest) velocity dispersion and, con- 
sequently, the least (most) concentrated density profile. 
The lightest star component is practically unaffected by 
the presence of the BH, while inside the BHIR the density 
cusp is particularly evident for the most massive stars. 

We note that the space density of each component 
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Table 2. Values for the components used in the multimass model. MS = ma in sequence 
stars ; G = giants; HB = horizontal branch stars; WD = white dwarfs (from lCote et all 
1 19951) . 



k 


mass range (Mq) 


M fe /M 


m fe (M Q ) 


(L/M) k 


o"fc/o" 


content 


1 


0.75-0.83 


4.9 x 10~ 2 


0.79 


10 


1 


MS, G, HB 


2 




3.9 x 10~ 2 


1.2 





0.81 


WD 


3 


0.65-0.75^ 


0.18 


0.70 


0.19 


1.1 


MS + WD 


1 


0.55-0.6~ 


9.3 x 10~ 2 


0.60 


6.5 x 10~ 2 


1.2 


MS 


5 


0.45-0.55^, 


0.24 


0.50 


2.3 x 10- 2 


1.3 


MS + WD 


6 


0.35-0.4~ 


0.16 


0.39 


1.0 x 10" 2 


1.6 


MS 


7 


0.25-0.35 


0.24 


0.29 


4.9 x 10~ 3 


5.7 


MS 



a Progenitors mass range that yield 1.2 Mq WD. 

^ Including a 0.7Mq WD population from progenitors with mass 1.5— 4 Mq. 
c Including a 0.5 Mq WD population from 0.83-1.5 Mq progenitors. 



must have the same asymptotic BW76 behaviour ~ x 7 , 
with 7 = 7/4 (i.e. ~ W 7 ^ 4 ), because it was assumed for 
simplicity that all the ft in equation (j 14|) have the same 
functional form. Nevertheless, how much deep one has to 
go towards the centre to find the x~ 7 ^ A cusp depends on 
Pk and, in turn, on the component stellar mass. This is 
because the density of the fcth component is a function of 
PkW oc a^ 2 [equation (|15[) ]; thus, the higher the velocity 
dispersion of the stars is, the deeper the extent of the 
region where the 'isothermal plateau' prevails, and the 
natter the density is at a fixed radius. 

With this in mind, in Fig. [6] the log slope 7 of the 
space mass density in the cusp region is plotted as a func- 
tion of the stellar mass. The slope is evaluated well inside 
the BHIR, namely at ~ xbh/5, i.e. where the enclosed 
mass in stars is of order 10 _2 M. . It can be noticed that 
sufficiently massive stars (m > 0.6 Mq) show a slope close 
to the linear relation 7 ~ 1.8m/(0.8 Mq), of the kind of 
that pred icted by the mult i com poncnt Fokker-Planck so- 
lutions of Bahcall fc Wola l| 19771 ) and then confirmed by 
iMurphv et all (|l99ll ). 

Given that low mass stars have a cusp region located 
deep within the core, t hey exhibit 7 < 1. 5 , in acc ordance 
to what was found in iBaumgardt et alj (|2004al ) . These 
authors also found a particularly low slope for intermedi- 
ate mass stars (figured)] dotted line), probably because in 
their simulations no particles r eached the x~ 7 ^ A re gime 
(see Fig. [S]and t he discussion in Frcitag ct al. 20063)- O n 
the other hand. iHopman fc Alexander 1 20061 ). including 
energy equipartition in the numerical integration of mul- 
ticomponent Fokker-Planck equations, reported power- 
law density behaviour with a high slope (~ 2) for massive 
objects (stellar BHs around the Galactic supermassive 
BH). 

For the mass range and mass function considered 
here, we found that the 7 - m relation is well fitted by 
the law 7 = [m/(l M ) - 0.29] 033 + 1. Although we can- 
not draw really general conclusions, it can be stated that 
our prescription to enforce energy equipartition generates 
results compatible with the findings of the various studies 
related to m ass segregation phenomena arou nd a massive 
BH (see also lAlexanderlll999l ; iMerritdliooj ). 

To give an indication on the chance of detecting the 



IMBH by means of observation of velocity dispersion pro- 
filefl let us suppose that the model generated here rep- 
resents a luminous and nearby cluster (with c ~ 1.2), 
located at a distance of 5 Kpc, with a central surface 
brightness Iq ~ 2 X 10 4 Lq pc -2 . If its core radius is a 
typical r c ~ 2 pc, the BHIR radius is ran = 0.094 pc. 
Moreover, let us assume that all the luminosity is given 
by the brightest component (that with fc = 1 in Table [2}, 
whose stars would have, at that distance, V ~ 16. The 
number of these stars inside the BHIR radius can be es- 
timated by 



Ni ~ /otitbh 



(£), 



mi 



(20) 



which gives Ni — 70, thus yelding relatively large Pois- 
son fluctuations that would make difficult to resolve the 
velocity cusp. 

This is confirmed by the example reported in 
Fig. [S^,, in which such measurements are done for a 
set of annuli around the cluster centre. The erro r bars 
are evaluated by the usual formulas (see, e.g., IJoned 
Il97d ; iMcNamara. Harrison fc "An derson 20Q3|) assuming 
an rms error on the single star velocity determination of 
~ 0.7 km s _1 (similar to, e.g., the typical accuracy level 
in the FL AMES-GIRAFFE V ery Large Telescope obser- 
vations bv lMilone et al.ll2006l . for stars with the same V) 
and treating the model a p profile as if it were the ob- 
served velocity dispersion profile with a typical a = 10 
km s" 1 . In spite of the relatively high Iq and low cluster 
distance, it is not straightforward to identify a cusp in 
velocity for such a non-collapsed cluster. Even if in a real 
observation the error bars could be smaller because of 
a larger number of stars actually observable, the detec- 
tion of a velocity cusp is certainly m ore difficult than for 
high-concentration clusters (see, e.g. , iGerssen et al II2OO2I ; 
iGersserJliobl ; iGebhardt et"afll2005j ). 

In general, from the grid of multimass models gen- 
erated, we found that ibh/ic ~ 0.07 ± 0.05, a ratio that 
will be used in the following section to give an estimate 
about Ni for the candidate clusters. 

4 Because of the assumed isotropy, line-of-sight velocity and 
proper motion dispersion are equivalent. 
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Figure 5. (a) Line-of-sight velocity dispersion profile and (b) 
projected surface density (normalized to the central value of 
the heavy remnants component) for the multimass model de- 
scribed in the text. Three components are shown: giant and 
turnoff stars (solid line), heavy dark remnants (dotted), light- 
est main sequence stars (dashed). The central logarithmic 
slope of £ are also reported for giants and heavy remnants. 
The vertical arrows indicate the BHIR radius. The error bars 
refer to simulated velocity dispersion measurements (see text). 
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Figure 6. Logarithmic slope of the space density of stars eval- 
uated for the model stellar components (dots) inside the cusp 
region (at x = 0.01z c ) and fitted by the law (m - 0.29) 33 + 1 
(solid line). The short-dashed line is the 7 /4 value predicted 
in the s ingle-mass context. The fit fro m the Baumgar dt et al.l 
J2004bl l results (dotted line) and the iBahcall fc Wolfl dl977T) 
relation (long-dashed) are also reported. 



4 BLACK HOLE MASS ESTIMATED FROM 
CLUSTER PROFILE 

Exploiting equation (|19[1 . we can now make a tentative 
estimate of the IMBH mass range in clusters, from their 
concentratiorjfl and from the slope of the central surface 
brightness profile. To t his pur poses, the recent and de- 
tailed analysis made bv lNGOtj of previous WFPC2/ HST 
observations of 38 galactic globular clusters, provides the 
necessary set of data. Moreover, we also considered the 
controversial case of the cluster Gl in M31, for which s 
is estimated by using t he recent ACS/ HST photomet- 
ric profile published in llvla et all (|2007l 'FI(that gives the 
most recently updated value for c, as well) and the total 
mass and luminousity taken bv lMevlan et al.l (|200ll ). 

From the whole list, 14 clusters were excluded 
(marked with ' x ' in the last column of Table [3} because 
they yield, from equation (|19[1 . a lower limit (M. t i) of the 
BH mass higher than the upper one {M.^)- All of core- 
collapsed clusters are in this subset, together with those 
having a steep central brightness profile (s > 0.2) in spite 
of being not very concentrated, like, e.g., NGC 6535 (that 
has c = 1.3 and s ~ 0.5). Of course, the presence of a 
IMBH cannot be completely ruled out for these clusters. 
It can only be stated that such a presence is incompat- 
ible with the stars phase-space distribution of the form 
given by equations ([T]) and (|14[l . Moreover, there are 18 
clusters (marked with 'o') having very low (or negative) 
central slope, for which M.,i < 100 Mq; conservatively, 
they were considered as clusters without IMBH. This be- 
cause the model is based on the assumption that M, is 
much greater than any single stellar mass, so that one 
can assume, in good approximation, the IMBH to be at 
rest. Therefore, as regards these clusters, the model can- 
not provide sufficiently stringent and reliable predictions; 
accurate individual parametric fit s would be necessary to 
this aim. Note, however, that the NG06 uncertainties on 
the measurements of s are rather large. This error is not 
taken into account in the present analysis. 

The relevant parameters of the remaining 7 can- 
didate clusters are written in Table 3] Their predicted 
IMBH mass range is plotted both vs. the total cluster V- 
band luminosity and vs. the observed velocity dispersion 
(Fig-0- It is apparent that M,/M has only a weak de- 
pendence on L while the BH mass tends to increase with 
<7 bs. To quantify these trends, weighted least -squares fits 
were computed (following the prescriptions of lPress et al.l 
1 19881 . sect. 15.2) neglecting, for simplicity, the uncertain- 
ties both in L and in a b s . The resulting fits are 

log (^f-) = (-1.2±1.6) + (-0.33±0.30) log f-^j (21) 

with xlt = 0.63 at a confidence level P(x 2 > Xflt) = 0.99, 
and 



Finally, note that the radius - measured in units of 
i c - at which the velocity dispersion is mostly affected 
by presence of the BH (where it follows a keplerian be- 
haviour) is ~ 0.02 for the luminous component (Fig. [5^0 . 
This factor is equal to abo ut 3M./M, in good agree ment 
with the relation found by I Baumgardt et al.l (|2005l ). 



5 According to what discussed in section [3] we can assume 
that c is equal to the observed standard concentration param- 
eter. 

6 s is evaluated using the central values of the surface bright- 
ness listed in their table 3 and a second order estimate of the 
first derivative. 



© ???? RAS, MNRAS 00Q.riTIT6l 



10 P. Miocchi 



log 



M. 

M 



(1.7± 1.4) + (1.2± 1.2) log (tJ^i) (22) 



with Xflt = 0.85 and P(x 2 > Xflt) = 0.97. The cluster 
luminosity was chosen as an independent variable because 
its measure is much more reliable and directly linked to 
observations than the d ynamically estimated total masses 
(|Prvor fc Mevlanlll993l ). 

Assuming a uniform global mass-to-light ratio 
equal to that averaged among the candidate clus- 
ters, ~ 3A(M/L) e , equation J2TJ gives M. ~ 
O.O9(M/M ) O 7±O - 3 M . This link is in agreement with 
the Mm oc M relation between super -massive BHs in 
galac tic nuclei and the bulge mass (|Magorrian et all 
1 19981 ). and between the mass of compact nuclei and that 
of the h ost spheroid in nucleated early-type and dwarf 
galaxies (|C6te et al.ll2006l ; IWehner fc Harrisll2006l ). 

As regards relation (|22[) , it is very different from 
the l aw M m ~ a^bs found b y iFerrarese fc MerrittJ 
(1200(f) for galactic bu lges (see also iGebhardt et al.ll200oT 
iTremaine et alj|2002j ). The reason of this discrepancy can 
be understood by the same argument IFerrarese fc MerrittJ 
presented to explain the relation they found: the M, 
- (Jobs link is a consequence of the fundamental rela- 
tion Mm oc M. In galaxies, this latter relation and the 
M ~ L 5//4 and L ~ cr^ ha laws lead indeed to M, ~ cr^ bs . In 
globular clusters, on the other hand, the observed trends 
M ~ L and L ~ a% 3 B l|Mevlan fc Heggid Il997h yield, 
through equation (|21|) . M, ~ fobsi which is, in fact, 
compatible with the relation (I22|l . 

However, since in old globular clusters the secular 
collisional relaxation has certainly had a substantial influ- 
ence on the dynamics within the core region, it is reason- 
able to expect only a weak correlation between M, and 
the presently observed core velocity dispersion. Moreover, 
as discussed in sect. 13.1.11 it is rather improbable that 
cr bs could provide a really reliable and accurate indica- 
tion on the velocity dispersion in the inner cusp region; 
in fact, with the exception of Gl, N\ is rather low (< 60) 
for the candidate clusters (see Table [4| . Conversely, the 
cluster total mass, in spite of the tidal erosion, should 
keep a tighter link with the environmental conditions at 
the first stage of the cluster life (at least for massive clus- 
ters), when IMBHs formation are thought to have oc- 
curred (e.g. lGiirkan et alj|200j ). 

To improve our capability of identifying the presence 
of IMBHs, it is important to search for other indicators 
that could be related to the gravitational influence of the 
compact object on the surrounding cluster stars. In this 
context, the presence of an extreme horizontal branch 
in the colour-magnitude diagram (CMD) may be a valid 
indicator, as we will see in the next Section. 

Not surprisingly, however, another indication in 
favour of the presence of IMBHs at the centre of our 
candidate clusters, is that their r c /r^ ratio - as directly 
evaluated from the values listed in the Harris catalogue, 
see Table [4]- is significantly higher t han t he maximum 
value (~ 0.05) suggested by iTrentil (2006) for clusters 
that do not harbour IMBHs, and, consequently, do not 
experience the enhanced core expansion caused by their 
presence. 

Note that three over five clusters that 




1 1.2 




log (L/L s ) 

Figure 7. Upper panel: BH mass of the candidate clusters 
plotted vs. the observed central velocity dispersion (in km 
s _1 ); lower panel: log of the ratio M,/M plotted against 
total cluster luminosity. The dashed lines are the weighted 
least-squares fits log(M./M) = -1.2 - 0.33 log(L/L ) and 
log(M./M ) = 1.7 + 1.21og(<7 ob s). 



iBaumgardt et all ((2005) indicate as possibly har- 
bouring IMBHs, are in common with our candidates 
set (namely: NGC 6093, 6266, 6388). We consider the 
remaining two clusters, NGC 5286 and 5694, as non- 
candidate because they have M m ,i > M % ,% (due to their 
relatively high concentration and slope); in this respect, 
notice that the core slope of the former cluster has been 
updated to an higher value in the fin al |NG06 | published 
list, compared with that reported by IBaumgardt et al.l . 
Furthermore, Gl, NGC 2808 and 6205 are not considered 
by these authors, because they identify as candidate 
clusters just those with a non-collapsed core and with a 
central slope between 0.2 and 0.3, as suggested by the 
results of their simulations. 

To conclude this section, it is worth noting th e simi- 
larity of the surface brightness profiles reported by |NG06l 
(plotted in their Fig. 7) with those illustrated here in 
Fig.s [TJj and 0). By inspecting those corresponding to 
our candidate clusters, the similarity is rather evident for 
M80 and M13. Indeed, these clusters show a clear transi- 
tion between a core a nd a cu sp region. Note that they are 
the only ones, in the |NG06 set, in which this feature is 
clearly observed. Thus, more detailed individual studies, 
also through parametric fits of the profiles based on the 
model discussed here, certainly deserve to be carried out 
for them. 



5 THE IMBH - EHB CONNECTION 

Let us assume that the loss of envelope mass is the dom- 
inant effect suffered by a giant during a close passage 
around the IMBH ('tidal stripping'), and that it be- 
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Table 3. Globular clusters in the lNGOrj set. L is the total luminosity in the V-band evaluated from My listed in the Harris 
catalogue, c = clus ter concentration (from the Harris catalogue; ':' denotes a core-collapsed cluster), s is the logarithmic 
slope ffrom lNGOrJ . with opposite sign). M #i i/M and M,^/M are the lower and up per limits of Mm /M as f rom equation 
1191 . M is the total cluster mass taken, if available, from the parametric estimates in lPrvor fc M cvlan ( 19931), while those 
marke d with a ':' ar e estimated by assuming a uniform global mass-to-light ratio = 2. EHBi marked clusters ar e those 
with a Harris (1996) HB type = 0. EHB2 clusters are those having log(T c g HB ) > 4.3 in the iRecio-Blanco et all ||2006T) 
list. Last column: • denotes a candidate cluster, 'x' a cluster having Mm 1 > M #j 2 and 'o' a cluster with M m ,i < 100 Mq. 
Not available data are indicated by .... 
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5272 (M3) 
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-2.41 


5.8 






• 


6254 (M10) 


4.92 


1.4 


-0.05 


-5.43 


-2.29 


5.4 









6266 (M62) 


5.61 


1.7 


0.13 


-3.34 


-2.63 


5.8 






• 


6284 


5.12 


2.5: 


0.55 


1.53 


-3.54 


5.4 






X 


6287 


4.88 


1.6 





-4.85 


-2.52 


5.2: 









6293 


5.04 


2.5: 


0.67 


2.92 


-3.54 


5.6 






X 


6341 (M92) 


5.21 


1.81 


0.01 


-4.73 


-2.76 


5.3 









6333 (M9) 


5.11 


1.15 





-4.85 


-2.00 


5.4: 









6352 


4.52 


1.1 


-0.02 


-5.08 


-1.95 


4.8: 









6388 


5.7 


1.7 


0.13 


-3.34 


-2.63 


6.2 






• 


6397 


4.58 


2.5: 
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7089 (M2) 
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-0.05 
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7099 (M30) 
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2.5: 
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6.31 


2.01 


0.06 


-4.15 


-2.99 


7.0 






• 



: and s comes from lMa et al. I 11200711 : L and M are taken from iMevlan et al] ll200lh . 



comes an EHB s tar because of i t s uncovered und erly- 
ing hotter layers (jRich et al ] |l997l : I Alexander! 120051 ). Let 
r tid = R,(M./m) 1/3 be the distance from the IMBH be- 
low which the tidal interaction with a star (with radius 
R, and mass m ) significantly perturb s its internal struc- 
ture (see, e.g., iFreitag fc Bend |2002| ; iBaumgardt et al.l 
2006). The order of magnitude of the rate of tidal strip- 
ping events can be estimated as S ~ noaoTvr^ id (l + 
2GM, /r t idO"o) 2 , where no and uq are, respectively, the 
typical stellar number density and velocity dispersion in 
the cusp region, and where the gravitational focusing is 



taking into account in the e stimate of the cross- section 
of the single stripping event l|Freitag et al. I l2006bh . 

For giants having R* ~ 100 R0 and m = mi, one 
has ryd — 2.5 x 10 -5 pc and taking, from the example dis- 
cussed in section I3XT1 n ~ 3A r i/47rr| H ~ 2 X 10 4 pc~ 3 , 
cr ~ 20 km s" 1 and M. = 10 3 Mq , one has S ~ 7 X 10~ 7 
yr~\ which means that of order of 100 stripping events 
can occur within a HB stars lifetime (~ 10 s yr), gener- 
ating a significant population of EHB stars. Incidentally, 
this population could have an integrated luminosity (10 2 
- 10 4 Lq, mainly in the UV band) compatible with the 
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Table 4. Candidate clusters set. Iq is the central surface brightness estimated from [iy inlfTarris ( f 996) using equation (5) in lDiorgovskil 
l|t993l l. r c is the core radius in pc, estimated from the angular value and the cluster distance listed in the Harris cat alogue. M»,i and 
M, 2 are in Mq and calculated using M from Table [3] cr b B is the central observed velocity dispersion taken from lPrvor &: Mevlanl 
l|l993l l or, marked with from Dubath, Mcvlan & Mavorj l|l997l ). M/L i s the mass-to-light ratio in solar unit. r c /r\^ is the ratio 
between the core and the half-mass radius, directly estimated from [Harris (1996). TVi is the estimate of the no. of giants inside xbh 
assuming ran = 0.07r c in equation (120 I t . The other symbols have the same meaning as in Table[3] 



NGC number 
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2808 


5.69 


3 x 10 4 


0.73 


110 


3.1 


X 
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13.7 ± 1.4 
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6093 (M80) 


5.22 
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0.44 
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10 3 
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13.4 ± 0.6: 


0.23 


11 


6205 (M13) 


5.41 


6.68 x 10 3 


1.7 


130 


2.4 
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10 3 


2.4 
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6.62 ±0.41 


0.52 


40 


6266 (M62) 


5.61 


2.54 x 10 4 


0.36 


290 


1.5 
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10 3 


1.5 




V 


14.3 ±0.4: 


0.15 


6 


6388 


5.7 


5.31 x 10 4 


0.35 


720 


3.7 


X 


10 3 


3.2 






18.9 ±0.8 


0.18 


13 


6715 (M54) 


5.94 


4.14 x 10 4 
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700 


3.2 
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10 3 


2.3 






14.2 ± 1.0 


0.22 


59 


Gl a 


6.31 


1.38 x 10 5 


0.78 


780 


1.1 


X 


10 4 


5.4 






25.1 ± 1.7 


0.12 


160 



and rjj comes from lMa et"al] {2007); M/L and <r obs are taken from lMevlan et"aL I ll200ll) 



far-UV excess showed by some - even metal-rich - clus - 
ters (e.g. NGC 6388, see lRich. Minniti fc Liebert]|l993l ). 
while, on the contrary, it would be insufficient if one con- 
sidered stripping events due to close encounters betwee n 
stars only, as already pointed out by iRich et al.1 1 19971 ). 
Thus, it is reasonable to expect that the existence of an 
EHB is favoured by the presence of an IMBH. 

To test this hypothesis, we rely on the classification 
made in the Harris catalogue (revised on February 2003) . 
This author re-a nalysed various set of CMD observations 
(mainly th ose by Piotto et al.l l2002). in the intent of up- 
dating the iDickensj (|l972F HB type, introducing a new 
type (= 0) to account for clusters having a prominent 
EHB. Thus, clusters with the HB type = are considered 
as having EHB (they are marked in the EHBi column of 
Table |3]) . Then, we study the probability distribution of 
finding a given number, Nehb, of clusters with an EHB, 
in a set as large as our candidate clusters set. However, 
for the sake of data homogeneity and because of the lack 
of direct observations of EHB stars, to the purpose of this 
analysis Gl is excluded from the candidates set. 

A large number of subsets made u p of 6 cl usters each, 
were extracted at random from the iNGOfj list. Then, 
-/Vehb was evaluated for each subset and a frequency 
distribution was computed (see the solid histogram in 
Fig- IEJ). The result is that the probability to get by chance 
a set of 6 clusters with JVehb ^ 3 (i.e. with A?ehb equal to 
or greater than that in our candidate clusters subset) is 
~ 0.063. This, from a statistical point of view, means that 
the hypothesis of the independence between the presence 
of an IMBH and the presence of an extended HB blue 
tail can be rejected at a confidence level > 90 percent. 

However, the Harris new HB type was assigned to 
clusters with an high relative proportion of blue subd- 
warfs population in comparison with that at the redder 
side of the RR Lyrae region of the HB (Harris, private 
communication), whileas the determination of the pres- 
ence of the EHB regardless of the other HB peculiar- 
ities would be more appropri ate to our purp o ses. T he 
various parameters defined bv lFusi Pecci et all (|l993l ) to 
characterize the HB morphology are unsuitable as well, 
because their sample contains only very few clusters of 



the ING06I list. For these reasons, we chose considering 
the maximum effect i ve tem perature, T e ffHB> evaluated by 
iRecio-Blanco et all (2006) for the HB stars of a clusters 
set that includes a significant part of the lNG06l sample 
(see their table 1). 

In this case, an EHB was considered as present when 
log(TeffHB) > 4.3 (i.e. TeffHB > 20,000 °K, see e.g 
Mochlcr 2003; iRosenberg. Recio-Blanco fc Garcfa-Marfnl 



20041 ). EHB clusters classified this way are signed 



in the ETB 2 column of Table [3] The visual inspec- 
tion of the IPiotto et al.| (|2002l ) CMDs - which the 
IRecio-Blanco et al.l ( 2006 ) evaluations are based on - 
confirms how in these clusters the EHB stars are clearly 
visible and reach downward a V magnitude as faint as 
the turnoff point. While a fraction of the clusters with 
such an EHB have an Harris HB type = 0, further clus- 
ters can be considered as owning EHBs. If one includes 
these clusters too, our candidates set has A^ehb = 4 and 
P{N E hb > 4) ~ 0.046 (see Fig. [8] dashed histogram), 
consequently the confidence level rises to ~ 95 percent. 

It is worth mentioning that the two candidate clus- 
ters (apart from Gl) considered not to posses an EHB, 
show significant peculiarities in their HB. In particular, 
NGC 6388 shows an over-luminous blue HB (IRich et al.l 
Il997l ; iRaimondo et all 120021 : ICatelan et ail 120061) and, 
simil arly to what observed in M54 1 Rosenberg et all 
2004|) , a blue HB exte nding anomalously beyond the the- 
oretical zero-age HB jBusso. Piotto fc Cassisill2004) . As 
regards Gl, we remind that there are some indirect indi- 
cations suggesting the presence of a non-n egligible pop- 
ulation of EBH stars l|Peterson et al.l [20031 ) and that, at 
least, such a presence cannot be e xcluded according to 
the most recent CMD observations (IRich et al.| [2005). 

It can be also noticed that in the lNGOd sample, there 
are 7 non-candidate clusters that possess EHB. This sug- 
gests that there might be other mechanisms - not related 
to the interaction with IMBHs - at work for the forma- 
tion of EHB stars. However, 4 of those clusters are highly 
concentrated (c > 2.4) or core-collapsed. This confirms 
that, whatever the formation mec hanism is, it is favoured 
by an high density en vironment (|Fusi Pecci et al] 1 19931 ; 
iBuonanno et al.ll 19971 ). Nevertheless, the presence of such 
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Figure 8. Probability distribution of finding TVehb clusters 
having an extended blue tail, in a subset of 6 clusters randomly 
chosen from the lNGoil set. The solid line histogram refers only 
to clusters with HB type = 0, while in the dashed one those 
having log(T c ff HB ) > 4.3 are also taken into account. Our 
candidate clusters subset has JVehb = 3 in the former case 
and = 4 in the latter. 



an environment at the present time is not a sufficient 
condition for the formation of an IMBH, which is closely 
linked to the initial central density instead. Gravother- 
mal oscillations and collapse may have changed rapidly 
the former dynamical state, thus erasing the link with 
the initial environment. A ll this could explain the ra ther 
weak correlation found by iRecio- Blanco et all i|2006l ) be- 
tween Teffjjg and the collisional parameter they denned 
on the basis of the present dynamical conditions. Fur- 
thermore, the velocity dispersion used in the definition 
could be not representative of the higher actual value in 
the inner (not resolved) cusp region. 

However, there are other alternative explanations 
for the presence of blue su b dwarfs, which are not based 
upon tidal stripping (see ICatelanl 120071 . sections 7.3 - 
7.4 and references therein). In this respect, notice that 
the helium self-enrichment - probably caused by multi- 
ple episodes of star formation - that in some clusters can 
explain t he existence of EHB populations (e.g., i n NGC 
2808, see lD'Antona et alJliooalPiotto et alj|2007l ). could 
be also related to the presence of an IMBH. In fact, the 
IMBH could originate subsequent bursts of star forma- 
tion, in a way similar - though on smaller scales - to what 
shou ld occur into the accretion disk of super-massive BH s 
(e.g. IShlosman fc Begelmadfl989l ; ICollin fc Zaimll 19991 ). 



6 CONCLUSIONS 

In this paper, the multimass isotropic an d spherical King 
model has been extended to include the iBahcall fc WolJ 
(| 19761 ) distribution function in the central region, so as 
to account for the presence of a central intermediate- 
mass black hole (IMBH). This self-consistent model gen- 
erates surface brightness profiles that show a power-law 
behaviour around the centre and a core-like profile out- 
wards. This latter is similar to a King profile with concen- 
tration c < 2, though slightly steeper in the core region. 
The following features are particularly relevant: 



c decreases monotonically with the IMBH mass M,; 

IMBHs with 10" 3 M < M m < 10 _2 M, with M the 
cluster mass, are compatible only with non-collapsed dis- 
tributions (c < 2); 

the core surface brightness profiles exhibit a logarith- 
mic slope s < 0.25 for reasonably low IMBH masses 
(Mm < 10~ 2 M). 

These features are in remarkable agreement with the 
results of recent collisional A-body simulations, which 
found that the IMBH causes a quick expansion of the a 
cluster core because of the enhanced collisional rate of 
exchange of kinetic energy both be tween stars of differ- 
ent masses l|Baumgardt et al"1l2005l ) and between single 
stars and hard binaries in 3-body encount ers occurring in 
the black hole vicinity dTrenti et al.ll200"7h . In particular, 
the slope of the profile of the core region in the simulated 
clusters attains values compatible with those found by 
our model. Moreover, the radius inside which the velocity 
dispersion profile is appreciably altered by the influence 
of the compact object (~ 3M,/M in u nits of the core ra- 
dius) coincides with the value found bv lBaumgardt et all 
We can also confirm that a full cuspy density profile (sim- 
ilar to, e.g., that of collapsed clusters) can be associated 
only to a black hole with a unrealistically high mass, even 
higher than that of the host cluster itself. 

A grid of models is then generated to derive pos- 
sible trends of the IMBH mass with the morphological 
parameters. It is found that, in general, 12s — 4.8 < 
log(M,/M) < —1.1c — 0.69. This mass range estimate is 
applied to a se t of 38 galactic globular clu sters recently 
re-analysed by iNovola fc Gebhardj (|2006h as well as to 
the case of Gl (in M31). It was found that seven clus- 
ters in this sample probably host an IMBH, namely: NGC 
2808, NGC 6388, M80, M13, M62, M54 and Gl. Among 
these, the M80 and M13 brightness profiles presented in 
INovola fc Gebhardtl show the typical appearance of the 
cuspy-core behaviour reproduced by the model. It is 
wort h noticing, moreo v er, th at among the five clusters 
that iBaumgardt et all (|2005h indicate as possibly har- 
bouring IMBHs, M80, M62 and NGC 6388 are in com- 
mon with our candidates set. 

The scaling relations 



Mm ~ 0.09 



M 



Mm ~ 50 



V km s 



M-. 



M,: 



(23) 



(24) 



are satisfied by our candidate clusters. Relation ([23|l 
is close to t he M, ~ M law found in larger 
scale systems (Magorrian et al.l Il99 1; ICote et all 120061: 
IWehner fc Harris! l2006h . None the less, relation (gljl is 
significa ntly different for the analogous relation found in 
galaxies llFerrarese fc Merrittll2000l ; lGebhardt et al.ll2000l : 
iTremaine et al.ll2002l ). but the discrepancy can be under- 
stood by noting that the scaling law between the total 
mass and the central velocity dispersion for galaxies is 
very different from that followed by globular clusters. 

Unfortunately, the uncertainties on the mass range 
estimates are in many cases rather large and, as a conse- 
quence, the presence of IMBHs on other clusters cannot 
be completely ruled out. They are mainly due to the un- 
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certainties in the central observed projected velocity dis- 
persion, and in the logarithmic slope measurements. For 
this reason, we plan to carry out in dividual studies on a 
repres entative cluster subset in the iNovola fc Gebhardtl 
( 2006) sample, to infer the IMBH mass directly from the 
best parametric fit of their brightness profiles. 

Nevertheless, this preliminary and 'collective' ap- 
proach, permits to achieve an intriguing result: the pres- 
ence of an extreme bl ue horizontal branch - determined 
according both to the iHarrisl (|l996T ) HB type index and 
to the maximum effective temperature of HB stars - is 
associated to the presence of the IMBH at a statistically 
significant level of confidence (> 90 percent). Such a cor- 
relation is not surprising when one takes into account 
the high stellar density in the inner cusp region and the 
relatively high ratio between M, and the single stellar 
mass, in calculating the rate of tidal stripping event in 
giants-lMBH close encounters. 

This firmly suggests that the mass of the IMBH 
could be one of the 'hidden' parameters that are being 
searched for to explain the strong variability of HB mor- 
phology in clusters with the same metallicity and similar 
CMDs (the so-called '2nd parameter' problem). For in- 
stance, the presence of a central IMBH in MI 3 and NGC 
6388 could explain the origin of their extreme HB stars, 
while their absence in the respective 'counterparts' M3 
and 47 Tuc is consistent with the the fact that, accord- 
ing to our analysis, IMBHs should not reside in these 
latter. 
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